High-Precision Measurement of the 19 Ne Half-Life and Implications for Right-Handed 

Weak Currents 
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We report a precise determination of the 19 Ne half-life to be T1/2 = 17.262 ± 0.007 s. This 
result disagrees with the most recent precision measurements and is important for placing bounds 
on predicted right-handed interactions that are absent in the current Standard Model. We are able 
to identify and disentangle two competing systematic effects that influence the accuracy of such 
measurements. Our findings prompt a reassessment of results from previous high-precision lifetime 
measurements that used similar equipment and methods. 

PACS numbers: 24.80. +y, 27.20. +n, 12.15.-y, 12.15.Hh, 12.60.-i, 29.40.Mc 



Precise measurements of decay rates and angular cor- 
relations in semi-leptonic processes are known to be ex- 
cellent probes for interactions that are predicted by ex- 
tensions of the Standard Model [lj. For example, the 
measured lifetime and electron asymmetry in neutron j3 
decay \2] are used to probe for right-handed currents and 
obtain a precise value of V u d, the up-down element of 
the Cabibbo-Kobayashi-Maskawa quark-mixing matrix, 
in a relatively simple system that is free of nuclear struc- 
ture effects. However, in spite of this compelling ad- 
vantage, precision neutron /3 decay experiments are chal- 
lenging. Current results from independent neutron decay 
measurements show large discrepancies that need to be 
addressed before conclusive interpretations can be made 
from the data 0. In this regard Nature offers a for- 
tuitous alternative in 19 Ne (3 decay, where the nuclear 
structure corrections are small and well understood 
Analogous to the neutron, the /3 decay of the isospin 
T = 1/2 19 Ne nucleus proceeds predominantly via a mir- 
ror 1/2+ — > 1/2+ transition to the analog ground state 
in 19 F. Furthermore, the vector and axial-vector form 
factors for the 19 Ne case are similar to neutron decay. 
These similarities make 19 Ne f3 decay studies an attrac- 
tive alternative to neutron decay experiments, without 
the experimental difficulties usually associated with the 
latter. It has been shown that the ft value (compara- 
tive half-life) for 19 Ne decay, together with the measured 
/3 asymmetry allows a determination of V u d H and 
constitutes one of the most sensitive probes for right- 
handed interactions in the charged weak current 6, 7]. 



Presently the dominant contribution to the total uncer- 
tainty in the measured ft value is from the 19 Ne half-life. 
Several measurements of the half-life have been reported 
that are in poor agreement with one another 

In this Letter we report a precise measurement of the 
19 Ne half-life, performed at the Triumf Isotope Separa- 
tor and Accelerator (Isac) facility in Vancouver, Canada. 
In this experiment, radioactive 19 Ne ions at ~ 37 keV 
were delivered to the 877 gamma-ray spectrometer, an 
array of 20 symmetrically placed Compton-suppressed 
high-purity germanium detectors. The inner volume 
of the array comprised of 20 similarly placed 1.6-mm- 
thick BC404 plastic scintillator detectors called Scep- 
tar (Scintillating Electron Positron Tagging Array) H, 
which were coupled to Hamamatsu H3695-10 photomulti- 
plier tubes (PMTs) and covered ~ 80% of the total solid 
angle. The ions were implanted on a ~ 1.3 cm wide, 
40 /im-thick mylar-backed aluminium tape at the center 
of the spectrometer, which allowed us to register data 
using tape cycles. In each cycle we took background 
data for the initial 2 seconds, collected the 19 Ne ions 
for ~ 1 second and reserved a counting time of 300 sec- 
onds (~ 20 half- lives). At the end of each cycle 1 second 
was allocated for the tape movement to move any exist- 
ing long-lived radioactivity away from the detectors, to a 
lead-shielded collection box placed downstream of the ar- 
ray center. The f3 signals from each Sceptar PMT were 
shaped by a fast preamplifier and input to a constant 
fraction discriminator (CFD) whose output pulses were 
sent to a logic OR fan-in/fan-out unit. Eight outputs 
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FIG. 1. (Color online) Summed TDC spectrum from all 
20 Sceptar detectors. The /3 — 7 coincidence, prescaled /3 
singles and the delayed prescaled /3 singles events are iden- 
tified. The CFDs were operated in "blocking" mode with a 
width of ~500 ns. This prevented them from re-triggering 
within 500 ns of the first pulse. The inset shows a close up 
of the afterpulsing events from the scintillator detectors and 
a fit to the data. 
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FIG. 2. (Color online) Extracted half-lives for a subset of 
the data (acquired at initial rates of ~1 kHz per detector, 
with 1 s bin-widths) corrected for different dead times. The 
upper panel shows a significant systematic effect that is more 
pronounced at lower dead times, where the fits to the data 
yielded worse x 2 ■ Normalized Monte Carlo simulation results 
are shown for comparison. The lower panel shows the same 
data set, corrected for afterpulsing using the results of our 
simulations. 



from the multi-channel fan-in/fan-out module were sub- 
jected to fixed non-extendible dead times in the range of 
4-24 /is. This time range was much longer than any other 
signal time in the electronics (500 ns from the CFDs), or 
the response of the pulses from the PMTs, thereby ensur- 
ing that pulse pile-up would have negligible effect. The 
counts from the dead-time-affected streams were regis- 
tered using a VME-based Struck SIS3801 multi-channel 
scaler (MCS) module. The time standard for the MCS 
data streams was provided by a 10 MHz ± 0.1 Hz temper- 
ature stabilized precision clock. Prescaled j3 singles (with 
a prescale factor of 255) and (3— 7 coincidence events were 
also recorded simultaneously in full list mode to obtain a 
measurement of the mirror (l/2 + — > 1/2+) branch. We 
defer discussion on this part of the analysis for a separate 
paper [§]. 

To test the ramifications of various systematic effects, 
we acquired the data at different initial rates and bin- 
widths, while the PMT bias voltages and CFD thresholds 
were kept unchanged. Two important systematic effects 
that influence the outcome of such high-precision life- 
time measurements are detector afterpulsing and PMT 
gain shifts. Afterpulsing depends on the type of scintil- 
lator and PMT used and leads to the detection of spu- 
rious counts so that the deduced half-life would appear 
to be lower than the true value. The low-amplitude af- 
tcrpulses mainly arise from long-lived phosphorescence 
in the scintillator material and the ionization of residual 
gas in the phototubes. On the other hand, PMT gain 
shift is a competing process that leads to event losses at 
high rates. These losses are caused by a fraction of the 
counts falling below the detection threshold and lead to 
an overestimation of the decay half-life. Below we de- 
scribe the methods used to minimize the contributions of 
both these systematic effects and the procedure used to 



identify and disentangle them. 

To minimize the effects of afterpulsing, we raised the 
CFD thresholds while simultaneously viewing the pulses 
on an oscilloscope. This procedure reduced the contri- 
bution of the afterpulses, without discriminating against 
real events. However, a time-to-digital converter (TDC) 
spectrum obtained from the event-by-event data stream, 
shown in Fig. [TJ revealed that a considerable number of 
the delayed pulses survived the raised CFD thresholds. 
The observed distribution of the remnant afterpulsing 
consisted of two exponentially decaying components with 
largely different decay constants. The areas under the 
prompt /3 singles peak and the delayed afterpulses showed 
the afterpulsing probability to be ~ 23%. The first step 
in the data analysis was to preselect the data by fitting 
the individual cycles from each experimental run to a 
function that consisted of an exponential plus a flat back- 
ground, using the maximum likelihood method described 
in Ref. The few cycles that yielded anomalously 

deviant half-life values were tagged and subsequently ex- 
cluded from further analysis. These comprised ~ 3% of 
the total data used for the final analysis. Data from the 
remaining 'good cycles' were then summed and fit using 
the same procedure. For dead time corrections to these 
data, the imposed dead times were determined using the 
source+pulser technique fllT | . The results from this anal- 
ysis are shown in the upper panel of Fig. [5J We obtain 
discordant half-life values for the same data that were 
corrected for different dead times. This is the first time 
that such an effect has been observed and reported in a 
scintillator-based, high-precision half-life measurement. 

To better understand this systematic effect, we per- 
formed Monte Carlo simulations of the events using the 
summed TDC spectrum as a model. In the simulation, 
the relative intensities and time distribution of the after- 
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FIG. 3. (Color online) Half-lives obtained by successive fits 
to data acquired at different rates. The bottom panel shows 
results from the summed data of four runs (representing all 
the low rate data) , while the other two panels show data from 
single runs. The bands in the high rate data show the time 
interval (after the beam was turned off) in which the rates 
would have dropped to ~ 450 Hz per detector. 



pulses that follow 19 Ne events were defined by the results 
of the fit shown in Fig. [TJ Simulated events that occurred 
outside a user-defined dead time window were registered 
separately. These data were binned, dead-time-corrected, 
and fit using the same method described above. The 
fits to the simulated data showed a similar underestima- 
tion of the extracted half-lives at the lower dead times. 
The general agreement between our model and the data, 
shown in the upper panel of Fig. [2j allowed us to make 
a first-order correction for the deviations. This was done 
using the following procedure: First, the simulated his- 
tograms were fit with an additional exponential compo- 
nent in the model, so that the modified fit function was 
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where At = — <; defines the width of the time bin. 
In the fits, the a 2 parameter was kept fixed at the input 
value for the 19 Ne half-life, while the other parameters 
were allowed to vary. The \ 2 probabilities of these mod- 
ified fits were found to be markedly better for the low 
dead-time data (ranging from 15% to 50%), with the fits 
converging to a value of a a — 8.81 ± 0.11 s at all dead 
times. The best-fit value for the a 3 parameter (which 
yields the relative intensities of the two components) de- 
pended on the initial rates input to the simulations and 
decreased as expected at larger dead times due to the 
presence of fewer afterpulses in the dead-time-corrected 
histograms. On fitting the experimental data similarly, 
with the 0,3 and 04 parameters kept fixed at the values de- 
termined by the simulations, we observe a restoration of 



the expected consistency (c.f. lower panel, Fig. [2]), which 
lends support to the validity of our simulations. Finally, 
additional simulations showed that the afterpulsing re- 
quired small corrections to the fixed dead times deter- 
mined using the prescription of Ref. [11]. These correc- 
tions (of the order < 1%) were incorporated accordingly. 

While the alignment of the extracted half-lives at dif- 
ferent dead times shown in Fig. [5] validates our under- 
standing of the afterpulsing, the accuracy of the absolute 
half-life is further influenced by other systematic effects 
that need to be taken into consideration. In the analysis 
described below, we investigate the other dominant sys- 
tematic effect due to PMT gain shifts. All fits described 
henceforth incorporate the corrections due to afterpuls- 
ing, using the function described in Eq. ([1]). We tested 
for rate dependent losses due to gain shifts by acquiring 
data at different initial trigger rates. As shown in Fig. |3l 
we observe significant losses at higher rates, which re- 
sult in the extraction of an erroneously large half-life. As 
an additional check for rate-dependent shifts, we succes- 
sively fit the histograms over a decreasing range of data 
points. In each fit, the initial time bin (that defined the 
first channel of the range over which the fit is performed) 
was increased in steps of 1-2 channels, and the half-life 
value was recorded. Fig. [3] shows the results from three 
scenarios, with data acquired at two extreme initial rates 
and an intermediate value of ~ 1 kHz per detector. We 
show only the results from the data that were subjected 
to the largest dead time of 24 f/s, with minimal afterpuls- 
ing correction. 

Although the data showed no obvious signs of rate- 
dependent shifts at initial rates of ~450 Hz, it was ap- 
parent that at the higher rates the PMT gains shifted 
significantly. Additionally, the data also showed that the 
phototubes needed a certain recovery time, depending on 
the intensity of the rates. This is obvious in the upper- 
panel of Fig. [3l where even after the detector rates had 
fallen to 450 Hz, the rate-independence shown in the bot- 
tom panel is not reproduced. Mindful of these effects we 
used only a subset of the total data from the experiment 
for the final analysis. These only included data that were 
acquired at initial rates of ^450 Hz and ~ 1 kHz per de- 
tector. For the data acquired at initial rates of ~1 kHz, 
we waited for a conservative time period of ~ two half- 
lives (34 seconds) after the beam was turned off to start 
fitting the data. With a relatively low background rate 
(~ 10 s -1 ), this approach ensured that the accuracy of 
our result was not compromised due to the recovery time 
of the PMTs mentioned above. 

Using the analysis described above we obtain the 19 Ne 
half-life to be T 1/2 = 17.262 ± 0.007 s. The individual re- 
sults from the runs are shown in the top panel of Fig. |4l 
where the central value in each run was obtained by tak- 
ing the mean of the half-life values obtained at the higher 
dead times of 12-24 (j,s (where the corrections due to af- 
terpulsing were the smallest). The contributions of the 
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FIG. 4. (Color online) Top panel: Extracted half-lives for dif- 
ferent runs at three different bin widths. The dotted line is the 
weighted mean, the lighter band shows the statistical uncer- 
tainty, and the darker band represents the final uncertainty in 
the half-life. Bottom panel: A comparision of the half-life val- 
ues obtained under different conditions to the weighted mean 
from all runs (with statistical uncertainties only). 



systematic and statistical uncertainties to our result are 
shown in Table HI The systematic uncertainty is dom- 
inated by gain shifts. We conservatively estimate the 
uncertainity due to gain shifts using the results obtained 
on varying the initial bin for the fits by ±12 seconds. 
Another systematic effect that could bias the result of 
this experiment is the potential long-lived diffusion of 
neon atoms from the tape. We assume this to be neg- 
ligible over the time intervals in this experiment. Our 
assumption is based on previous measurements, where 
no apparent neon diffusion was observed using the same 
tape [l2j . The bottom panel in Fig. @] shows that our 
measurement is not significantly affected by other sources 
of systematic uncertainties. 

Our result for the 19 Ne half-life disagrees with 
the most recently published high-precision result, 
17.219±0.017s [13], by approximately 2.5 standard de- 
viations. Using the present value of the half- life, the Qec 
value of the decay (liL th e analog branch Q and the elec- 
tron capture fraction 15j, |16[ , we obtain the corrected Tt 



TABLE I. Uncertainty budget for the half-life measurement. 



Source 


Uncertainty(s) 


Afterpulsing corrections T 


0.003 


Gain shifts 


0.006 


Dead time corrections 


0.0005 


Statistics 


0.002 


Total 


0.007 



Obtained from our Monte Carlo simulations. 
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FIG. 5. (Color online) In the above, the bands represent the 
±68% confidence intervals. Gv = GpVudgv, where gv is 
the vector form factor. In the low momentum transfer limit, 
gv(q 2 -»■ 0) = C v = 1. 



value of this mixed transition using the equation 

Tt = .f v t{l + 5' R ){l + 5 v NS -5 v c ) 
K 



G 2 y M|(l + A^)(l + ^) 



(2) 



= 1721.3 ± 1.2 s. Here G v = G F V ud C v is the effec- 
tive vector coupling constant, and we follow the defi- 
nitions and notation described in Refs. 0, 0]. The re- 
sulting Tt value is nearly a factor of two more precise 
than the previously adopted value Together with 
the vector coupling constant, this value can be used to 
predict the beta asymmetry for 19 Ne decay assuming 
only V — A interactions and the conserved vector cur- 
rent hypothesis. We find the Standard Model beta asym- 
metry parameter to be = -0.0416 ± 0.0007, which 
can be compared to the experimentally measured value, 
A® xp = -0.0391 ± 0.0014 Q. This comparison is shown 
in terms of the mixing ratio and the vector coupling con- 
stant in Fig. [5j where the constraints set by our half-life 
measurement can be seen. The data show a hint of dis- 
agreement which invites further investigation. A remea- 
surement of the beta asymmetry parameter (which dates 
to 1975) is welcomed at this point. 

In conclusion, our measurement provides useful data 
to complement ongoing precision neutron decay experi- 
ments and is particularly relevant for searches of right- 
handed interactions in the charged weak current. We ob- 
serve a combination of systematic effects that need to be 
considered for an accurate result. We speculate whether 
similar effects, particularly detector afterpulsing, could 
have gone unnoticed in other high-precision half-life mea- 
surements that used plastic scintillators. 
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